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We propose that dramatic changes in the variances and covariance of protons and antiprotons can 
result if baryons approach chemical equilibrium in nuclear collisions at RHIC. To explore how equi- 
libration alters these fluctuations, we formulate both equilibrium and nonequilibrium hadrochemical 
descriptions of baryon evolution. Contributions to fluctuations from impact parameter averaging 
and finite acceptance in nuclear collisions are numerically simulated. 
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Event-by-event fluctuations of the particle yields in rel- 
ativistic nuclear collisions can be sensitive to the degree of 
chemical equilibration |jj and to critical fluctuations due 
to the QCD phase transition |2j|j| . We study for the first 
time the degree to which chemical equilibration can affect 
the fluctuations of baryons at RHIC and LHC. Baryons 
and antibaryons are likely produced in such abundance 
that event-by-event yields can be measured Q. We pro- 
pose that chemical equilibration can appreciably change 
the variances and covariance of protons and antiprotons 
compared to thermal and participant-nucleon model ex- 
pectations. 

Fluctuations of the baryon-antibaryon system can vary 
depending on the degree to which the chemical reactions 
NN <-» mesons reach equilibrium In the absence 
of chemical equilibrium, the number of antibaryons and 
baryons are independently conserved. Fluctuations of the 
particle number are then Poissonian, 



In chemical equilibrium, we find that the covariance can 
be negative, 



V = ^(-/Vj — (N)) 2 ps TV, thermal 



(1) 



where Ni is the baryon number in the i th event and 
N = (N) is the event average. HIJING and similar 
participant nucleon models exhibit essentially the same 
behavior for nucleus- nucleus collisions In contrast, 
chemical equilibrium relates fluctuations of the baryons 
to those of the antibaryons, reducing the relative variance 
compared to (Q). We find ||: 



V ps N 2 (N + Al- 



chemical 



(2) 



where N is the average number of antibaryons. We re- 
cover (|j) in the limit N 3> N, since baryon conservation 
fixes the number of particles in that limit. 

If chemical equilibrium is achieved in RHIC colli- 
sions, the variance (^) can fall to half the thermal and 
participant-nucleon model level, because the numbers of 
baryons and antibaryons are expected to be compara- 
ble at midrapidity M. We propose that a further sig- 
nal of equilibration can be obtained from the baryon- 
antibaryon covariance, 



C ps -NN(N + Al- 



chemical 
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C = X>,- (N)Wi- (AT)). 



(3) 



On the other hand, this quantity essentially vanishes for 
collisions of large nuclei in thermal equilibrium, or when 
N > AT. ' 

The anti-correlation (^) seems somewhat surprising - 
positive correlations between baryons and antibaryons 
are seen in e + e~ experiments [§] and incorporated, e.g., 
in string fragmentation models JtJ. Such correlations 
are required by baryon number conservation for those 
comparatively small systems. In contrast, the anti- 
correlation (||) and the reduction of the variance (||) result 
from fluctuations of the net baryon number itself. In a 
subvolume of a large system, fluctuations can increase the 
net baryon number Nb if either a particle enters or an 
antiparticle exits the subvolume. Correspondingly, fluc- 
tuations of Nb simultaneously increase A?" and reduce AT, 
leading to an anti-correlation. This effect is only possi- 
ble if chemical reactions render the individual numbers 
A" and N indefinite. 

To determine the degree of chemical equilibration from 
fluctuations in experiments, one must account for the 
following additional sources of fluctuations, all of which 
tend to increase V and drive C toward positive values. 
Near local equilibrium, thermal and volume fluctuations 
produce particle number fluctuations. In experiments, 
further fluctuations result from centrality selection and 
detector acceptance effects HH. To illustrate how vol- 
ume and thermal fluctuations can change the variance 
of particle yields, we begin by extending the thermody- 
namic formulation of ref. Jpf. We then apply the Monte 
Carlo event generator of Qto compute the experimental 
contribution. 

Our results imply that novel correlations can be ob- 
served if chemical equilibrium is obtained. However, it is 
not obvious that baryons will achieve chemical equilib- 
rium or that equilibrium correlations will survive chem- 
ical and thermal freezeout. We develop a simple hadro- 
chemical model to study the approach toward chemical 
equilibrium. We estimate the corrections to (|l|) and (||) 
that arise if the system is only partially equilibrated. We 
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then briefly discuss how elastic scattering following chem- 
ical freezeout can alter C and V 

To illustrate how fluctuations alter antibaryon produc- 
tion in local thermal and chemical equilibrium, we em- 
ploy the idealized but standard Bjorken hydrodynamic 
framework. We suppose that the energy and baryon num- 
ber are deposited near midrapidity by pre-equilibrium 
processes that locally establish an initial temperature T, 
entropy density s and net baryon density ps- These 
quantities are further assumed to vary only with proper 
time r. Entropy and baryon number conservation then 
imply that the rapidity densities for entropy, S = dS/dy, 
and net baryon number Nb — dAfs/dy are constant. 
The rapidity density of all hadrons is also nearly con- 
stant, because iVtot oc S for a system dominated by light 
hadron species with masses <C T. Mean rapidity den- 
sities of individual species, such as antibaryons N and 
baryons N vary with r. For this Bjorken scenario, it 
is useful to define an effective volume V w At that in- 
creases from a formation time tq to freezeout at Tp, such 
that V — S/s oc N tot / n tot . The transverse area A is ini- 
tially determined by the overlap of the colliding nuclei. 

In local thermal and chemical equilibrium, the vari- 
ances and covariances of the state variables specify the 
fluctuations of all bulk thermodynamic quantities. For 
our collision scenario, the appropriate state variables are 
T, V and Nb- The fluctuations of a quantity X are 
characterized by the variance o\ = (AX 2 ) for AX 
X — (X). If we take the total number of hadrons Not to 
be proportional to the entropy, then 



r 2 y/V 2 



(AN, 
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As in we take the fluctuations of T and Nb to satisfy: 
4 = T 2 C-\ a% = TdN B /dp B (6) 



for a system with a heat capacity C v 
For an ideal hadron gas, 



see e.g. ref. 
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T 2 ay /12V' 2 , a B =N + N, 



(7) 



neglecting small corrections from Fermi and Bose statis- 
tics. To completely specify our ensemble, we assume that 
fluctuations of these state variables are statistically inde- 
pendent, so that (AVAT) = (AN B AT) = (AN B AV) = 
0. Observe that while this formulation gives the follow- 
ing the flavor of a thermodynamic analysis, ours is not a 
static, global equilibrium scenario. 

Volume and thermal fluctuations cause the number of 
baryons and antibaryons to fluctuate B . In the absence 
of chemical equilibrium, we write: 



AN = NAV/V + AN V 



(8) 



with a similar equation for antibaryons. The contribution 
AN VT for constant T and V is ([!]). Thermal fluctuations 
do not change the conserved particle numbers, while vol- 
ume fluctuations add to ([I]). The variance is 



Vth = = N 2 a 2 v /V 2 + N^N(\ + N/N tot ). (9) 

We also compute the baryon-antibaryon covariancc: 

C th ee (AN AN) = NNa 2 v /V 2 m NN/N tot . (10) 

This result implies that given N particles, the probability 
of finding an antiparticle is N/N to t- The correction to p) 
and ([!]) from volume fluctuations are of order 4%, since 
N/Ntot ~ 0.04 for central Au+Au at RHIC. 

In chemical equilibrium, the number of particles and 
antiparticles obey: 



AN = 



and 
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AT+ (S) ^ (12) 



where (8N/dT) NB = (dN/dT) NB and (dN/dN B ) T - 
(dN / ON b)t = 1. In chemical equilibrium, thermal fluc- 
tuations can change the baryon population by making 
pairs, so that 



dN 



For an ideal gas, 



dN f dN \ dN B /dT 



dT \dfiB J dNs/nB 



AN = N 



AV 2eNN AT NAN f 



V N+N T N+N 



and 



AN = N 



AV 2eNN AT NAN F 
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N+N T N+N 



(13) 



(14) 



(15) 



where e = En /NT w m/T + 3/2 is the energy per an- 
tibaryon per unit temperature. We then obtain: 



V ch = N 2 ^ + 4e 2 



N 2 



NN 



N + N J T 2 N + N' 



For N S> N we obtain (g). The covariancc is: 

NN \ 2 al NN 



ch 



V 2 



NN-^r + 4e 2 



N + N J T 2 N + N 



(16) 



(17) 



The first terms in ( |m are the volume contribu- 
tions encountered earlier. These terms are unchanged 
by chemical equilibrium. The second term accounts for 
the pairwise production of baryons and antibaryons by 
thermal fluctuations. The third term describes the net 
baryon number fluctuations discussed earlier, see eq. (Q). 

At RHIC, N w N implies that both the baryon vari- 
ance and the baryon-antibaryon covariance markedly dif- 
fer from the nonequilibrium results. The sum of the vol- 
ume and thermal terms is roughly 1 + e 2 /12 w 5.5 for 
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chemical freezeout at T = 160 MeV and m = 938 MeV. 
The contribution from the first two terms to either 
V ch /N 2 or C ch /NN is then - 0.5% for iV tot - 10 3 as 
expected in Au+Au collisions. Baryon density fluctu- 
ations contribute ~ +1.3% to the variance and —1.3% 
to the covariance for A w A ss 40, yielding the totals 

V ch /A 2 ~ 1.8% and C ch /NN 0.8%. Observe that 

the chemical nonequilibrium results (|[ |h}) give values 
~ 2.5% and +0.1% for the relative variance and covari- 
ance; we expect participant nucleon models to yield sim- 
ilar values. We will see that that these small percentage- 
differences amount to quite large changes in the magni- 
tudes of these quantities. 

We comment that particle ratios, such as A/A to t, are 
often measured in place of absolute yields. Variances 
of ratios do not contain the contribution from volume 
fluctuations found e.g., in @ and (16). However, ex- 
perimenters must carefully construct ratios from A and 
iV to t from data measured over the full kinematic ranges 
of these particles in order to fully cancel the volume fluc- 
tuations and to suppress new fluctuations due to system- 
atic shifts in the spectra. The distinction between ratios 
and yields is minor, however, since we find that volume 
fluctuations are the smallest of the contributions treated 
here. 

To account for additional sources of fluctuations in 
heavy ion collisions, we incorporate these general results 
into the Monte Carlo event generator for correlated sig- 
nals developed in Q). There, we assume that the rapidity 
densities A, A and A to t can be described by a thermal 
ensemble at each impact parameter b, with fluctuations 
occurring about the well defined mean values obtained 
below. We then generate events, choosing the values of 
N l and Ni for the i th event using (§), Q or (§]]), 
depending on whether the system is in thermal equilib- 
rium, chemical equilibrium, or in between. To simulate 
the centrality dependence of ion collisions, we distribute 
events in b according to a wounded nucleon model with 
realistic nuclear density profiles. 

We compute the mean rapidity densities of protons 
and antiprotons for Au+Au at RHIC using the wounded 
nucleon model relations, 

A « rjjV(b)/A/(0), A m nM(b)/JV{0), (18) 

where Af(b) is the number of participant nucleons. We 
take values of the rapidity density at zero impact param- 
eter n « n « 40 for both baryons and antibaryons. Such 
values are consistent with the range of event generator 
predictions: HIJING, HIJING/BB (its successor) and 
RQMD report rapidity densities of 60, 20 and 20 respec- 
tively. [Note that we will use (|l^) to provide initial con- 
ditions when we consider partial equilibrium effects, see 
eq. (22).] We then compute the average charged particle 
multiplicity for each event assuming a Gaussian distribu- 
tion with an average value N tot (b) = N tot (0)N(b)/N(0) 
and a standard deviation <r tot = \/Atot consistent with 
thermal equilibrium. The scale A to t(0) = 2100 is deter- 
mined by the initial production regardless of event class, 



in accord with entropy and energy conservation; the par- 
ticular value is taken from a HIJING simulation in the 
STAR acceptance. 

The negative covariance that signals chemical equi- 
librium is not destroyed by volume, thermal, finite- 
acceptance, or impact-parameter fluctuations. In fig. 1, 
we show the proton variance and proton-antiproton co- 
variance computed from 10 6 events in the STAR accep- 
tance assuming that chemical equilibrium is complete. 
Such an event sample can be accumulated in roughly 
twelve days of STAR running and can provide a statisti- 
cally significant determination of C and V as functions of 
the multiplicity. Chemical equilibrium fluctuations are 
computed using ([11, 12)- Our results shown by the solid 
curves are strikingly different from thermal equilibrium 
expectations - the long-dashed curves - obtained from 
eq. (||) at the same density. Incorporated in our simula- 
tion is an important non-thermal source of fluctuations 
that results from impact parameter averaging, as dis- 
cussed in 0]. These "background" fluctuations, shown 
as the short-dashed curves, result from the binning of 
the variances and covariance as functions of multiplicity. 
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FIG. 1. Proton variance (top) and proton-antiproton co- 
variance (bottom) for a chemical equilibrium state in RHIC 
Au+Au. The signature of baryon-density fluctuations - 
the negative correlation - remains striking in the covari- 
ance, despite competition from thermal, volume, and im- 
pact-parameter fluctuations. 

Having found that chemical equilibrium for baryons 
can have observable consequences, we now ask whether 
these particles approach chemical equilibrium in the first 
place. To be concrete, we suppose that baryons, an- 
tibaryons and other hadrons form at a proper time tq 
from the hadronization of a quark gluon plasma. These 
baryons can then scatter and annihilate, while reactions 
such as puj — > NN produce new baryons. Chemical equi- 
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librium is obtained when the annihilation rate precisely 
matches the creation rate. 

To describe the approach to equilibrium, we extend the 
kinetic model of |^| to include meson-meson production in 
the following schematic fashion. We take the antibaryon 
and baryon density to satisfy the approximate kinetic 
equation: 



d 

dr' 



1 



+ - ) n = -(a a V Te i) (nn - n cq n cq ) 



together with baryon number conservation. Again we 
assume a Bjorken-like expansion. Detailed balance for 
TV TV # mesons fixes the product of the chemical equilib- 
rium densities n cq n cq in terms of the meson densities. 
Note that this relaxation-time approximation strictly 
holds only when the initial baryon, antibaryon and meson 
densities are sufficiently close to local chemical equilib- 
rium values. 

In terms of the rapidity densities TV and TV, we write 



rdN/dr fa -7 {NN - TV cq TV cq }. 



(20) 



Baryon conservation implies N — N = TVg. The coeffi- 
cient 7 = (cr a w r ei)^4 _1 depends on the collision frequency. 
In H , we used the measured energy-dependent annihila- 
tion cross section to compute (a a v TC i) ps 4.4 fm 2 . We 
take iVeq to be a r-independent parameter to be varied 
over a range of plausible values. This assumption is rea- 
sonable for N, N <C iVtotj since the change in the meson 
populations due to baryon annihilation is then negligible. 



the initial rapidity density TVo satisfies Nq < TV eq , i.e., if 
the ratio 



IV„, 



N 



N eq +N -N B 



is greater than zero. The solution is then: 

n-^ + (n 

2 + I cq 2 1 + 9S 



(21) 



(22) 



(19) 

where 



s = (To/T F yV N ««- N °\ 



(23) 

with chemical freezeout occurring at rp . If 9 <C 1 then 
annihilation reduces the antiproton-to-proton ratio as in 
The evolution of N toward chemical equilibrium is 
shown in figs. 2a and 2c for TV cq > No and TV cq < No, 
respectively. 

To estimate fluctuations in the nonequilbrium state, 
we write the following ansatz: 



V = SV th + (1 - S)V, 



ch, (24) 

where Vth and V c h respectively satisfy (§) and (|l^). To 
motivate this ansatz, we consider the variance as a func- 
tion of t and linearize for N near TV c h, as is standard 
0. We find V - V ch + 2TV cq (TV(r) - N cq ). Differenti- 
ating and using (|2p|), we see that rdV/dr pa 2kV, where 
k ~ d(\ogN)/d(\ogT) - -7. The ansatz @ has the 
correct long- and short-scale time dependence to linear 
order, satisfies the appropriate boundary conditions for 
t near zero and infinity and has plausible behavior for 
intermediate r. The ansatz 



C — SCth 
is similarly plausible. 



(1 - S)C 



ch 5 



(25) 
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FIG. 2. The rapid approach to chemical equilibrium of 
the proton rapidity density (top) and variance (bottom) as 
functions of the freezeout time tf for N cq > TVo, (a) and (b), 
and for N cq < Nq. Hadron formation occurs at to = 8 fm. 

To solve (^o|), observe that N(t) will increase toward 
N cq depending on whether the right hand side of ( |2p| ) is 
positive or negative. Growing solutions are obtained if 




FIG. 3. Proton variance (top) and covariance (bottom) 
for a partial chemical equilibrium state in RHIC Au+Au com- 
pared to complete chemical (dot-dashed) and thermal equi- 
librium (long dashed) expectations for the same N and TV. 
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To appreciate the speed with which equilibration takes 
place, we show the calculated approach to chemical equi- 
librium of iV and V in fig. 2. We assume that hadroniza- 
tion proceeds from a thcrmalized quark gluon plasma 
and, correspondingly, take To = 8 fm. The initial values 
Nq « Nq are taken from (|l8|) for 6 = 0. The chemical 
freezeout time is varied up to a value Tp ~ 14 fm, a plau- 
sible value for central Au+Au at RHIC. For this figure 
we take ad hoc values for N eq « iVeq to illustrate the 
approach to equilibrium from above and below. Equi- 
libration is rapid because the annihilation cross section 
and baryon rapidity densities are sufficiently large that 
the ratio of the collision time to the expansion time scale 
satisfies T coU /r cxp ~ (27A cq ) _1 < 1. 

The variance and covariance computed for the STAR 
acceptance using 10 6 events simulated according to our 
partial equilibrium solution are shown in fig. 3. Rapidity 
densities are computed using (|l] - ^3|) for initial rapidity 
densities from (jig). In these calculations, we take the 
ad hoc value N cq = 50 for central collisions and assume 
N eq scales with the number of participants. The chemi- 
cal freezeout time is computed for each impact parameter 
using tp = R/v s + t p , where R = (A/it) 1 / 2 is the geo- 
metric transverse radius, v s = l/v3 is the sound speed 



and 



1 fm is a limiting value for peripheral collisions 



set by the proton size. Using these values, fluctuations 
are simulated in accord with |||). We again take 
the formation time for hadrons to be To = 8 fm inde- 
pendent of b. Note that for values r ~ 5 fm or smaller, 
chemical equilibration is practically complete for impact 
parameters smaller than 10 fm. 

We now ask if such fluctuations can survive freeze- 
out. If chemical freezeout occurs before thermal freeze- 
out, then baryons can suffer elastic collisions with other 
hadrons after chemical fluctuations are no longer pos- 
sible. Such scattering can restore C and V to thermal 
equilibrium values. To estimate the size of this effect, we 
assume that those baryons and antibaryons that do not 
scatter elastically contribute to fluctuations at the level 
of V c h and C c h, while baryons that scatter once or more 
are instantly thermalized. We then write Jl2| , 



v ~ rv ch + (i - v)v. 



th 



(26) 



where the survival probability for elastic scattering is 



V 



- \ dt av YC Yn to t{r+vt,t) 



(27) 



for n tot the hadron density, a the elastic scattering cross 
section and v re \ the relative velocity. Pions dominate 
the late-time dynamics, so that a k, a n N ~ 20 mb. 
To estimate V, we assume that after chemical freeze- 
out the meson and baryon distributions are described 
by spherically symmetric Gaussian profiles for both po- 
sition and velocity. Averaging over these distributions, 
we find 0.8 < V < 0.9 for central collisions, assuming 
TttotCrp) ~ 0.1 fm~ 3 and 1 < v s t f /R < 2. 

These calculations suggest that scattering after chemi- 
cal freezeout increases the chemical equilibrium variances 



and covariance by less than 20%. Scattering at that level 
is not likely to wipe out the large chemical fluctuations in 
figs. 1 and 3. Furthermore, UrQMD calculations suggest 
that chemical freezeout for baryons at RHIC may occur 
much later than we have assumed pl| . In that case, we 
expect rescattering to be completely negligible. 

In summary, we have found that baryon density fluc- 
tuations lead to striking correlations of protons and an- 
tiprotons if RHIC collisions approach chemical equilib- 
rium. Measurements of the variance and covariance of 
these species can reveal these correlations, even if equi- 
libration is incomplete and chemical freezeout precedes 
thermal freezeout. We have not examined how the p 
and p measurements are altered by resonance decays or 
hyperon contributions (see [fl3[ for a discussion of res- 
onances in another context). In addition, we have ne- 
glected mean field and coulomb effects, which modify 
correlations over times 3> Tp. Such effects are impor- 
tant when the momentum dependence of correlations is 
resolved in an HBT-like analysis, see e.g. ]l4| and refs. 
therein. Our effect contributes pre-freezeout correlations, 
often neglected in HBT discussions. 
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